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Abstract 


If massive neutrinos are Dirac particles, the proposed PTOLEMY experiment will 
hopefully be able to discover cosmic neutrino background via ve +°H —> ?He+e~ witha 
capture rate of Tp ~ 4 yr~'. Recently, it has been pointed out that right-handed com- 
ponents of Dirac neutrinos could also be copiously produced in the early Universe and 
become an extra thermal or nonthermal ingredient of cosmic relic neutrinos, enhancing 
the capture rate to Tp © 5.1 ye? or Ip + 6.1 yr !. In this work, we investigate the 
possibility to distinguish between thermal and nonthermal spectra of cosmic relic neu- 
trinos by measuring the annual modulation of the capture rate. For neutrino masses 
of 0.1 eV, we have found the amplitude of annual modulation in the standard case is 
M x 0.05%, which will be increased to 0.1% and 0.15% in the presence of additional 
thermal and nonthermal right-handed neutrinos, respectively. The future detection of 
such a modulation will be helpful in understanding the Majorana or Dirac nature of 


massive neutrinos. 
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1 Introduction 


The discovery and a subsequent precision measurement of cosmic microwave background 
(CMB) have provided us with precious information on the evolution of our Universe, and 
helped us to establish the standard model of cosmology [1]. In contrast, only an indirect 
evidence for the existence of cosmic neutrino background (CvB) has been found in the 
observations of CMB, primordial abundances of light nuclear elements and cosmic large-scale 
structures, where neutrinos have played a crucially important role [2] B, 4 [5]. The direct 
detection of CvB in the future will be another great triumph of the standard cosmology, and 
also offer a new possibility to probe the intrinsic properties of massive neutrinos. 

According to the standard cosmology, the temperature of CvB at present is extremely 
low, namely, T, ~ 1.95 K, and the corresponding average momentum is (p,,) ~ 5x 10~* eV ØB] 
[4 B]. Consequently, it will be an extremely difficult task to detect CvB, although the 
number density for each neutrino or antineutrino species can be as large as 7, œ~ 56 cm7°. 
As was pointed out by Weinberg in 1962 [6], cosmic relic neutrinos can be captured on 
beta-decaying nuclei, such as v, + 3H — 3He + e~, which can take place without an energy 
threshold for the initial-state neutrinos. In Ref. [7], it was noticed that the signal of the 
capture of CvB will be a peak in the final-state electron spectrum, which is located at a 
distance of 2m,, away from the end point of beta-decay spectrum of 7H — *He+ 7, + e. 
Therefore, in order to observe this tiny separation, we need an energy resolution even smaller 
than the absolute neutrino mass m,,, which is now known to be at or below the sub-eV level. 

Recently, a lot of attention has been paid to the detection of CvB through the capture on 
beta-decaying nuclei. See, e.g., Refs. [8] and [9], for reviews on the other possible approaches. 
In Ref. [10], a detailed calculation of the capture rate for v, +°H — *He + e~ was first 
performed, and the impact of neutrino mass hierarchies and flavor mixing on the capture 
rate was investigated in Ref. [11]. In addition, similar processes for the eV- and keV-mass 
sterile neutrinos or for different nuclear isotopes have also been discussed in the literature [12] 
PIJ. More importantly, the PTOLEMY experiment has been 
proposed [22], in which a 100 gram of surface-deposition tritium source will be implemented 
and the energy resolution for the final-state electrons can reach A = 0.15 eV. Given this 
nominal experimental setup, Long et al. and Lisanti et al. have found that the capture rate 
is Tp ~ 4 yr! if massive neutrinos are Dirac particles, while that for Majorana neutrinos 
is twice larger, i.e., Ty œ% 8 yr™! (24). Although the capture rate is small, it is already 
interesting to have a realistic experiment that will be able to directly prove the existence 
of CvB. Moreover, the difference between the capture rates I and Fy will enable us to 
discriminate Dirac and Majorana nature of massive neutrinos. Such a discrimination is 
extremely important for exploring the origin of neutrino masses. 

Motivated by the proposed PTOLEMY experiment, the authors of Refs. and have 
examined how the copious production and subsequent decoupling of right-handed neutrinos 
in the early Universe can affect the capture rate of cosmic relic neutrinos, assuming massive 
neutrinos to be Dirac particles. It has been shown that the capture rate can be increased at 
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most to be Tp © 5.1 yr7! or Tp & 6.1 yr~', if right-handed Dirac neutrinos are produced 


thermally or nonthermally [26], respectively. In this framework, an immediate question 
arises: is it possible to distinguish between thermal and nonthermal energy spectra of relic 
right-handed neutrinos? The main purpose of this short note is to answer this question, and 
show that a PTOLEMY-like experiment with a total statistics of 10° neutrino events may 
do such a challenging job. 

The remaining part of this work is structured as follows. In Section 2, we briefly review 
possible ways to generate right-handed neutrinos in the early Universe and the impact of 
thermal and nonthermal spectra of relic right-handed neutrinos on the detection of CvB. 
Following the general idea of Ref. [27], where the gravitational focusing effects by the Sun 
are considered, we demonstrate in Section 3 that the annual modulation of the capture rate 
can be used to discriminate the velocity distributions of cosmic relic neutrinos. Finally, we 


give some concluding remarks in Section 4. 


2 Cosmic Relic Neutrinos 


If the standard model of elementary particles (SM) is extended with three right-handed 
neutrinos, which are singlets under the SM gauge symmetries, neutrinos acquire tiny Dirac 
masses in the same way as charged leptons and quarks do. More explicitly, the gauge- 


invariant Lagrangian can be written as 


L=Lyoy t+ Y Tiluan- |S) Yap lor Arth], (1) 
a a,b 

with a and @ running over e, u and 7, where Lg), stands for the SM Lagrangian, laL 
and H = io,H* are lepton and Higgs doublets, Y, is the Dirac neutrino Yukawa coupling 
matrix. After the spontaneous gauge symmetry breaking, the neutrino mass matrix is given 
by M, = Y,(H) with (H) ~ 174 GeV being the vacuum expectation value of the Higgs field. 
In order to obtain sub-eV neutrino masses O(M,,) ~ 0.1 eV, we have to require extremely 
small Yukawa couplings O(Y,) ~ 107!?, which exaggerates the hierarchy problem of fermion 
Yukawa couplings, when compared to the top-quark Yukawa coupling y, ~ 1. On the other 
hand, the lepton-number conservation or a global U(1) symmetry should be imposed on the 
Lagarangian in Eq. to forbid a Majorana mass term for right-handed neutrino singlets. 
In the following discussions, however, we focus on the cosmological implications of massive 
Dirac neutrinos. A satisfactory solution to these two theoretical problems is beyond the 
scope of this work. 

Within the SM and standard cosmology, the production rate of right-handed neutrinos 
in the early Universe is always far below the Hubble expansion rate at any time [25]. 
Therefore, it is quite reasonable to ignore right-handed neutrinos vp and left-handed an- 
tineutrinos Vp for a = e, u,T in the evolution of our Universe. Nowadays, we have only 
left-handed neutrinos v, and right-handed antineutrinos 7,p as a cosmic background, for 
which the average number density is 7, ~ 56 cm~’. When new physics scenarios are con- 
sidered, it is possible to produce both v,p and Y, abundantly and reconcile them with 
cosmological observations, in particular the effective number of neutrino species N gẹ from 
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the latest data of Big Bang Nucleosynthesis (BBN) and CMB. Generally speaking, there 
exist two typical scenarios. 


e Thermal Right-handed Dirac Neutrinos—As is shown in Ref. [25], if primordial mag- 
netic fields are generated in the electroweak phase transition [29], even a rather small 
magnetic dipole moment of massive Dirac neutrinos u, = 3 x 107 up (with upg being 
the Bohr magneton and m, ~ 0.1 eV) leads to a thermal production of vp via signif- 
icant Var, > Var conversions [80]. Depending on the evolution of primordial magnetic 
fields, the decoupling temperature of v g can be well above the temperature of QCD 
phase transition Tocp © 200 MeV. Hence, a proper decoupling temperature Tp can 
be chosen to evade the cosmological bound N ş = 3.144933 at the 95% confidence 
level [31]. In this scenario, vąę and Va, will survive today as extra thermal relics, 
whose distribution functions in the CyB frame are isotropic and can be parametrized 
as follows 

1 


exp(p,/TR) + 1’ 


fm) = (2) 
where Tf is the present temperature of vag and Yar. For the decoupling temperature 
Tr © 200 MeV, the cosmological bound on AN, = Nog — 3.046 < 0.53 is saturated 
and we have TR ~ 1.28 K. The average number density of v,, or Tr turns out to 
be 7, (var) ~ 16 cm~’, which should be compared with 7,(v,,,) ~ 56 cm~? in the 
standard case. Thus, in the thermal (TH) scenario, only Tf < 1.28 K is allowed by 


cosmological observations. 


e Nonthermal Right-handed Dirac Neutrinos—In Ref. [26], it has been argued that the 
coupling of right-handed neutrinos to the inflaton field could give rise to a nonthermal 
production of v,p and Ya in the reheating epoch after inflation. In this nonther- 
mal (NT) scenario, both vg and F,r, have never been in thermal equilibrium, the 
distribution function in the CvB frame can be written as 


n, p, SEE; 
fyr(p,) = (3) 


0, p, >e, 


where £% denotes the Fermi energy at present. Here 0 < 7 < 1 is the fraction of 
occupied states below the Fermi energy, and 7 = 1 corresponds to the case of a degen- 
erate spectrum. Without specifying any dynamics of particle production from inflaton 
decays, the authors of Ref. [26] have assumed that the Fermi energy at production is 
Ep = ETR, where Tp is the temperature of thermal bath and € is a free parameter. 
The cosmological bound AN, < 0.53 can be translated into a restrictive constraint 
on this free parameter € < 3 for a degenerate spectrum with 7 = 1] If this upper 
bound is saturated, today we have the Fermi energy €} ~ 2.75 K, and the average 


'To be consistent with the TH scenario, we use the same upper bound AN.g < 0.53 on extra neutrino 
species, while AN. < 0.7 was adopted in Ref. [26]. Therefore, a slight difference between the number 
density T, (Vag) ~ 29 cm~® in the present work and 7, (Vap) ~ 36 cm~® in Ref. [26] should be noted. 
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Normalized Velocity Distribution 
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Figure 1: The normalized velocity distribution gy(v,) = v2g(v,)m3/(27?7,,) for thermal and 
nonthermal relic right-handed Dirac neutrinos in the CvB frame, where the temperature 
Tg = 1.28 K is taken for the former case and the Fermi energy £% = 2.75 K for the latter. 
Note that the absolute neutrino mass m,, = 0.1 eV is assumed for illustration and the speed 
of light c = 3.0 x 10° km s~ is used. 


number density of var or Tar is T, (Var) = 29 cm, which should be compared with 
Ti (VaL) ~ 56 cm~® in the standard case. For a nondegenerate spectrum with 7 < 1, 
the number density will be slightly smaller. 


Notice that the cosmological bound on AN. < 0.53 will be further improved in the future. 
For instance, if the upper bound on AN, is reduced by a factor of two, the allowed number 
density of right-handed neutrinos 7,(V,p) will be multiplied by 278/4 = 0.59 for both TH 
and NT cases. 

Throughout this work, we shall take the absolute scale of neutrino masses to be m, ~ 
0.1 eV, which is marginally compatible with the cosmological bound X) = m, + m, + 
Ma < 0.23 eV on the sum of neutrino masses f] Hence, cosmic relic neutrinos are nowadays 
non-relativistic particles. For later convenience, we write the distribution function of non- 
relativistic relic neutrinos as 
Gp, 3 dt 


f(e) r)? 7 (anya) (4) 


where p, ~% m,v, holds as an excellent approximation and g(v,) is the velocity distribution 
function in the CvB frame. 


?It is worthwhile to mention that a more stringent bound © < 0.12 eV or < 0.13 eV at the 95% confidence 
level has been derived by combing the Planck CMB and Baryon Acoustic Oscillation data with the Lyman- 
alpha clustering data or with the SDSS-DR7 sample of Luminous Red Galaxies [33]. Once such a tight 
bound is confirmed, the PTOLEMY experiment will have to greatly improve its energy resolution to well 
separate the signal from a huge background. In addition, neutrino masses below m,, œ% 0.1 eV imply larger 
neutrino velocities, reducing the annual modulation, as we will discuss in the next section. 


The integration of the right-hand side of Eq. over neutrino velocities gives us the neu- 
trino number density 7,,. In order to compare between TH and NT scenarios, we introduce 
the normalized velocity distribution as gy(v,) = v2g(v,)m3/(27°7,,) and show it for right- 
handed neutrinos in Fig. 1, where one can observe the main difference between these two 
types of spectra. The NT spectrum concentrates more on the velocity region of 200 km s~! 


to 700 km s~! in the CvB frame, while the TH spectrum extends to a much higher velocity. 


3 Annual Modulation 


In the standard case, it has been proposed in Ref. [27] that the gravitational focusing of CvB 
by the Sun will modify the distribution function of v,;, at the Earth and cause an annual 
modulation of the capture rate of CvB in the PTOLEMY experiment. We shall apply this 
interesting idea to the relic right-handed Dirac neutrinos vp and examine if it is possible to 
distinguish between TH and NT spectra, arising from different production mechanisms. In 
this section, we first summarize the capture rates of cosmic relic neutrinos for a PTOLEMY- 
like experiment. Then, following Ref. [27], we introduce the measure of annual modulation 
of the capture rate for standard CyB and compare it with those for STD+TH and STD+NT 
with relic right-handed Dirac neutrinos] 


3.1 Capture Rates 


The capture rate of a polarized neutrino state on the tritium target was first correctly 
calculated in Ref. [23]. See Ref. [24] for a similar calculation, and also some discussions on 
measuring neutrino anisotropy on the polarized tritium target. Considering an unpolarized 
tritium target and a neutrino mass eigenstate v; with a spin s, being either +1/2 or —1/2, 
one can find the product of the cross section and the neutrino velocity Vy, 


o,(s,)v,, = A(s,)|U4|?F ’ (5) 


where F ~ 3.8x 1074 cm? and A(s,) = 1—2s„v, , and U,, for i = 1, 2,3 are the three elements 
in the first row of the lepton flavor mixing matrix. Since A(+1/2) = 0 and A(—1/2) 2 hold 
in the relativistic limit v, >, only the left-helical neutrinos can be captured. However, in 
the non-relativistic limit v, — 0, we have A(+1/2) ~ A(—1/2) © 1, so both left- and right- 
helical neutrino states contribute to the total capture rate. Summing over the contributions 
from all three neutrino mass eigenstates, we obtain the total capture rate [23] 


3 d3 
Po =N fae (pafa Ped + OE afaa Ped] > (6) 


where N, stands for the number of tritium nuclei, Jas (p, ) and Fox (p, ) denote the distri- 


bution functions of left-helical and right-helical neutrino states at the Earth, respectively. 


3To avoid any confusion, we emphasize that TH and NT are referring only to right-handed neutrinos, 
while STD+TH and STD+NT to both left- and right-handed neutrinos. 


In the standard scenario of massive Dirac neutrinos, the production of right-handed neu- 
trinos Vp in the early Universe is so inefficient that all the neutrino states are left-handed in 
chirality. In the epoch of neutrino decoupling, neutrinos are relativistic particles, so almost 
all of them are left-helical states. As the helicities are conserved during the expansion of our 
Universe, we have the number density of left-helical neutrino states ni ~ 56 cm7? nowa- 
days, where all the neutrino mass eigenstates are assumed to be equally populated. Mean- 
while, the number density of right-helical neutrino states is vanishingly small, i.e., Ty x0. 
Following the same arguments, we arrive at ns ~ 56 cm? and ny x 0 for antineutrinos. 
Since g;v, 1s independent of neutrino velocities in the leading-order approximation for non- 
relativistic neutrinos (i.e., v, & 1), the integration over distribution functions just gives the 
neutrino number density. With the help of the unitarity condition of lepton flavor mixing 
matrix |U,4|? + |U.|? + |U.3/? = 1, the capture rate of v, + °H — He + e7 at PTOLEMY is 
found to be Tp ~ 4 yr™' B3]. 

In the nonstandard scenario of either TH or NT production mechanism, the right-handed 
neutrinos Vag and left-handed antineutrinos 7,;, can be populated, although their number 
densities are well constrained by cosmological observations. In assumption of equal number 


3 


densities for neutrino mass eigenstates, we obtain at most Mia = ni ~ 16 cm™ and 


Mi = ne ~ 29 cm~? for TH and NT cases, respectively, depending on the decoupling 
temperature Tp in the former case while the initial Fermi energy €p and the occupation 
fraction 7 in the latter. Neglecting the velocity dependence and inserting these number 
densities into Eq. (6), one can get the total capture rate Tp © 5.1 yr7! and Fp & 6.1 yr7! 


for STD+TH and STD+NT scenarios, as already mentioned before. 


3.2 Focusing Effects 


To discriminate the velocity distributions of relic neutrinos, as shown in Fig.|1| now we should 
take into account all possible dependence on neutrino velocities when calculating the capture 
rates. First, the spin factor A(s,,) in Eq. actually depends on neutrino velocities, namely, 
A(s,) = 1 +v, for left-helical states and A(s,) = 1 — v, for right-helical states |] Although 
such corrections to the total rate are insignificant, we include them for completeness. 
Second, the impact of the Sun’s gravity on the dark matter distribution has already been 
investigated in Refs. [B4 [36]. As has been pointed out in Ref. [27], even if the relic 
neutrinos are unbound to the Milky Way, the gravity of the Sun will modify the velocity 
distribution of CvB at the Earth, where the detectors are located. In connection with 
gravitational effects on neutrinos, it is worthwhile to notice that massive neutrinos are subject 
to gravitational clustering in the dark-matter halo of the galactic scale, which has been 
studied carefully by following the nonlinear evolution of CvB in Refs. B38]. For m, = 
0.1 eV the relic neutrino overdensity n,,/m, — 1 is found to be more than 10% in Ref. [38], 
while for m, = 0.15 eV a value about 40% is reached in Ref. [37]. Thus, it may be reasonable 
to expect an enhancement of 30% in the local neutrino density. On the other hand, some 


4If massive neutrinos are Majorana particles, the distribution functions of left-helical and right-helical 
neutrino states are identical and thus the velocity dependence in A(s,,) is completely cancelled out [27]. 
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Figure 2: The normalized capture rate ',(t)/T'§'? (tmin) of relic neutrinos in the PTOLEMY 
experiment has been shown for the standard (STD) scenario (blue solid curve), the nonstan- 


min) 


dard scenario of thermal (TH) right-handed Dirac neutrinos (orange dotted curve) with 
Tg = 1.28 K, and the case of nonthermal (NT) right-handed Dirac neutrinos (red dashed 
curve) with e% = 2.75 K., where P3!°(¢ 


case. 


min) denotes the capture rate at tmin in the standard 


deviations from the Fermi-Dirac distribution have also been observed [37 [38]. As shown in 
Ref. [38], the proportion of neutrinos with peculiar velocities below 100 km s~! is smaller by a 
factor of 1.33 than that for the Fermi-Dirac distribution. The impact on higher velocities will 
be much smaller. Though we ignore these effects in our calculations, a relative error about 
30% should be kept in mind when we calculate the capture rate and annual modulation. 

In order to derive the distribution function of relic neutrinos relative to the Earth, which 
will be denoted as gẹ(v,), we have to know the velocity of the Earth Vẹ relative to the Sun. 
Working in the heliocentric ecliptic coordinate system, for which the ecliptic plane is spanned 
by two unit vectors ê and ê, one can find Vz, = Vz [ê; cosw(t — tye) + êa sin w(t — t,,)] 
with V} = 29.79 km s“! and w = 2r yr7', where t,e March 20 is the time of the vernal 
equinox [36} [39]. The velocity of the Sun relative to the CvB frame is approximately given by 
Vo = Vo(—0.0695, —0.662, 0.747) with Va = 369 km s“!. The velocity Vo can be identified 
from the dipole structure of CMB |40] if we approximately identify the CMB frame with 
CvB frame, which is quite reasonable as no significant clustering of neutrinos takes place in 
the epoch of CMB formation [27]. Finally, due to the Sun’s gravity, the neutrino velocity v, 


at the Earth corresponds to the velocity v,, at an infinitely far distance to the Sun 


KaD va V, + Vvo(GyuMo/lr|)f = Dsl Vy ` r)v 


= - í 
v,,(V_) v2, + GyMo/lr] — vlv,- Î) ’ ( ) 


where v = |v,,| and f = r/|r| with r being the displacement vector of the Earth, Gy 
Newton’s gravitational constant, and Mọ the solar mass. According to the Liouville theorem, 
the distribution function is constant along any trajectory in the phase space. Thus, we have 


9e(V.) =9[Vot+VoolVetv)] ; (8) 


for the local distribution function of relic neutrinos at the Earth. As is shown in Ref. [27], 
the minimal neutrino density in the standard case is reached around tmin = tey — 8 days 
(namely, March 12) when the Earth is most upwind, while the maximum is achieved half a 
year later, i.e., on September 11. 

Now we can insert the distribution function g,(v,,) for relic neutrinos in Eq. together 
with Eq. into the total capture rate in Eq. (6), and derive the rates at an arbitrary time 
during the year both in the standard scenario and in two nonstandard scenarios with right- 
handed neutrinos. The final results of [p(t)/T3"? (tin) are illustrated in Fig. |2| where all 
the total capture rates [5 (t) have been normalized to the minimal rate [3'? (t 


min) 
min) on March 
12 in the standard case. Since O{Vy, is almost velocity-independent, the normalized capture 
rates directly reflect the local neutrino number densities along the orbit of the Earth. 

Two main differences between standard and nonstandard cases can be observed. First, the 
relic right-handed neutrinos in STD+TH and STD+NT scenarios contribute to the capture 
rate by a sizable amount of 28% and 52%, respectively. Second, the difference between 
maximum and minimum increases remarkably in the two nonstandard cases, where right- 
handed neutrinos concentrate more on the region of low velocities and thus they are affected 
greatly by the Sun’s gravity. To be explicit, the average neutrino velocity in the STD case 
is approximately 1589 km s~! for m,, = 0.1 eV, while the average velocities of right-handed 
neutrinos for TH and NT cases are around 1043 km s~' and 533 km s~, respectively. So the 
right-handed neutrinos in the STD+TH and STD+NT scenarios will increase the fraction 


of low-velocity neutrinos, when compared to the STD case. 


3.3 Detailed Comparison 


To make a detailed comparison, we follow Ref. and define the annual modulation as 


_ Tp(t) —Lp (tem) 
es as OES CHW 9) 


where tnin © March 12 stands for the time of minimal capture rate. In Fig. |3| the annual 


modulation has been calculated for both standard and nonstandard scenarios. On the left 
panel, the modulation for the standard CvB is plotted as a magenta dotted curve, from which 
one can read the maximum is only about 0.05% (See also the same curve zoomed in on the 
right panel). In the TH case, the thick and thin dashed curves correspond to the modulations 
for the temperature of T} = 0.3 K and TR = 1.28 K, respectively. Similar results for the NT 
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Figure 3: The annual modulation M(t) of the CvB capture rate in the PTOLEMY exper- 
iment, where the modulation for the standard case (STD) is represented by the magenta 
dotted curve, while those for TH and NT cases by the dashed and solid curves, respectively. 
On the left panel, the modulation is shown only for right-handed neutrinos, and two differ- 
ent values TR = 0.3 K and TQ = 1.28 K are chosen for the TH case while e% = 1 K and 
ey = 2.75 K for the NT case. On the right panel, the contributions from both left-handed 
and right-handed neutrinos to the capture rates are combined. 


case of ef = 1 K and eb = 2.75 K are denoted by thick and thin solid curves. On the right 
panel, the combined contributions from both left-handed and right-handed neutrinos to the 
modulation have been depicted. Some comments are in order: 


e One can observe from the left panel of Fig. [3]that the peak modulation can be as large 
as 3.5% in the TH case for T} = 0.3 K and in the NT case for e£% = 1 K, which should 
be compared to M < 0.05% in the standard case. For the other model parameters 
Tg = 1.28 K and £% = 2.75 K, for which the number density of right-handed neutrinos 
saturates the cosmological upper bound on AN, the maximal modulations in both 
TH and NT will be reduced by a factor of twenty. This observation demonstrates that 
the impact of the Sun’s gravity on the local density of relic neutrinos is more important 


for slower neutrinos. 


e In reality, the modulation should receive the contributions from both left-handed neu- 
trinos (i.e., the standard CvB) and right-handed neutrinos, and one cannot turn off the 
former contribution. Therefore, only the curves labelled by STD+TH and STD+NT 
on the right panel correspond to realistic situations. On the left panel, the modula- 
tions solely from right-handed neutrinos are given to illustrate their main differences 
from the result in the standard case. Although the amplitude of annual modulation 
due to right-handed neutrinos is larger for a lower temperature or Fermi energy, their 
number density will be smaller, leading to a reduced capture rate. This is why the 
peak modulations for different model parameters (namely, TR for TH and £f for NT) 
on the right panel are now comparable. 
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Figure 4: The amplitude of the modulation M is shown with respect to the ratio of right- 
handed neutrino number density Np to the photon number density N_. The results for only 
right-handed neutrinos in two nonstandard scenarios are given on the left panel, while the 
combined contributions from both left-handed and right-handed neutrinos are shown on the 
right panel where the corresponding capture rate solely from right-handed neutrinos is given 
on the top axis. 


Finally, we consider if it is possible to distinguish between STD+TH and STD+NT scenarios, 
assuming the number density of right-handed neutrinos to be the same. This consideration 
does make sense, since the total capture rates will be distinguishable if the number densities 
of right-handed neutrinos are different. Hence, in Fig. |4| we show the maximal modulation 
for a given number density Np of right-handed neutrinos, which has been normalized to the 
present-day number density of CMB photons N, ~ 411 cm~°. In either TH or NT scenario, 
once a value of T$ or £} is assumed, one can immediately calculate the corresponding number 
density Np and the capture rate PR, which are actually displayed on the bottom and top 
axes. Some important observations related to Fig. |4|should be emphasized: 


1. For an extremely small number density Np, the contribution of right-handed neutri- 
nos to the total capture rate is negligible, as is shown on the right panel of Fig. 
Therefore, the modulation for STD+NT (blue solid curve) and that for STD+TH (red 
dashed curve) are approaching M ~ 0.05% from right, which is the limiting value in 
the standard case. To the right of the vertical dashed line around Ng/ N, x 107, 
it is evident that the modulations for STD+TH and STD+NT scenarios are clearly 
distinguishable. 


2. It should be noted that the slowly-moving neutrino wind does not imply a large mod- 
ulation for extremely low neutrino velocities. The main reason can be understood as 
follows. In the heliocentric ecliptic coordinate system, the wind of CvB is moving 
towards the solar system in an inclined angle about 0.067 or 10.8°, which indicates 
that the focus point through the solar gravitational lens does not exactly match the 
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maximum point of modulation. Because of the dispersion of neutrino velocities, the 
focus of neutrino wind is broadened to cover the maximum point of modulation. It is 
straightforward to verify that this can be achieved if the neutrino velocity increases to 
Vo x 0.067 ~ 70 km s71. Therefore, if neutrino velocities are too low, the focus region 
is far from the orbit of the Earth, implying a small modulation. 


3. For TH and NT right-handed neutrinos, their number densities Np can be calculated by 
using the distribution functions in Eq. and Eq. (3), respectively. Thus, the average 
temperature TR or Fermi energy £% of right-handed neutrinos in the TH or NT case 
is proportional to NĚ 3 implying that a larger value of Np corresponds to more faster 
neutrinos. On the left panel of Fig. |4| one may notice that there exists a maximum 
of the modulation for varying number densities Np in either NT or TH scenario. The 
reason is the same as above. The average velocity at the maximum point is 67 km s7! 

for the NT case, while 73 km s~! for the TH case, which is very close to 70 km s7! 
that we have found before. After these maximum points, the modulation decreases due 
to the fact that the neutrino wind is less deflected for increasing velocities. The two 
vertical lines on the right part arise from the cosmological bound on AN, constraining 
the number density of right-handed neutrinos from above. On the right panel, however, 
those maxima are no longer visible, which can be understood by the relatively small 
number density of relic right-handed neutrinos at the maximum compared to that of 
left-handed neutrinos. Note that new maxima indeed arise from the superposition of 
the contributions from both left-handed and right-handed neutrinos and are located 
at different number densities. 


It is now clear that the maximal modulation in the nonstandard STD+TH and STD+NT 
scenarios ranges from 0.1% to 0.15%, which is twice or three times larger than M œ~ 0.05% 
in the standard case. In order to observe such an effect, we need about 10° events, which can 
be accumulated in a PTOLEMY-like experiment with a tritium target of two ton running 
for ten years. In accumulating the data, we should also notice that the incoming direction 
of the neutrino wind is changing, leading to a drift of the maximum point of modulation at 
about 0.3° per 20 years [39]. 


4 Concluding Remarks 


The recently proposed PTOLEMY experiment will be able to detect relic neutrinos from the 
Big Bang via ve, +°H — *He+e7 by using a 100 gram of surface-deposition tritium source and 
an excellent energy resolution of 0.15 eV [22]. Assuming that massive neutrinos are Dirac 
particles, the authors of Ref. [25] and Ref. [26] have found that the capture rate of cosmic 
relic neutrinos can be enhanced from Tp ¥ 4 yr7! in the standard case to be Ty œ% 5.1 yr! 
and Ty © 6.1 yr~', if right-handed neutrinos can be produced thermally and nonthermally 
in the early Universe. In this short note, following the idea of annual modulation induced by 
the Sun’s gravity in Ref. [27], we have pointed out that it is possible to distinguish between 


thermal and nonthermal spectra of cosmic relic neutrinos. 
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In the optimistic case, where the cosmological bound on the extra neutrino species 
AN. < 0.53 is saturated, the modulation of the CvB capture rate reaches its maximum 
M = 0.15% or M & 0.1% on September 11 for the STD+NT and STD+TH scenario, 
respectively. This is to be compared with the maximal modulation M œ 0.05% in the stan- 
dard case for the absolute neutrino mass scale m, = 0.1 eV. The choice of this absolute 
mass is marginally compatible with the cosmological bound on the sum of three neutrino 
masses X < 0.23 eV. For a larger neutrino mass, the amplitude of modulation will be 
increased and the gravitational clustering effects will enhance the local number density of 
cosmic relic neutrinos. In this case, however, neutrinos will be a mixture of unbound and 
bound components, implying an even smaller difference between STD+NT and STD+TH 
spectra. As shown in Ref. [27], the bound neutrino wind has a different focal point around 
the Sun, implying a different phase of annual modulation. The velocity distribution of either 
NT or TH right-handed neutrinos will be affected significantly by gravitational clustering 
in the dark-matter halo of our galaxy and the difference between them will be diminished. 
Therefore, a relatively large neutrino mass is not always helpful. 

As a rough estimate, we have to accumulate more than 10° events in a PTOLEMY-like 
experiment to see an annual modulation of 0.1% at the 20 level [Z7]. Such a detection might 
not be impossible in the future when the PTOLEMY experiment is successful in discovering 
cosmic relic neutrinos and then scaled up to a 10 kg target mass or more. Furthermore, the 
general idea presented in this note is to discriminate thermal and nonthermal energy spectra 
of cosmic relic particles via the annual modulation of detection rates. It can be applied to 


e eV-mass sterile neutrinos — They are able to explain a possible extra radiation during 
the BBN and CMB era, and can be implemented to relax the tension among cosmo- 
logical parameters determined from data sets [41]. Such light sterile neutrinos can be 
produced in the early Universe via a mixing with active neutrinos, or exotic interactions 
with dark matter. 


e keV-mass sterile neutrinos — They are one of the most promising candidates for warm 
dark matter [42]. In the early Universe, there are several different ways to produce 
keV-mass sterile neutrinos, including a mixing with active neutrinos in the presence of 


a large lepton number asymmetry and decays from a singlet scalar. 


For these sterile neutrinos, a small mixing with active neutrinos guarantees a similar capture 
process and an annual modulation at the PTOLEMY experiment. However, the masses of 
sterile neutrinos are much larger than those of ordinary neutrinos, so it may be important 
to investigate their evolution in the gravitational potential and figure out their velocity 
distributions in the present solar system. As both the production mechanism and thermal 
history of such particles can be very different from ordinary neutrinos, we may even have a 
distinct phase in the annual modulation. We leave these interesting topics for a future work. 
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